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Abstract 

In theories with a gauge-mediated mechanism of supersymmetry breaking the 
gravitino is likely to be the lightest superparticle and, hence, a candidate for dark 
matter. We show that the decay of the next-to-lightest superparticle into a gravitino 
can yield a non-thermal population of gravitinos which behave as a hot dark matter 
component. Together with the warm component, which is provided by the popula- 
tion of gravitinos of thermal origin, they can give rise to viable schemes of mixed 
dark matter. This realization has some specific and testable features both in par- 
ticle physics and astrophysics. We outline under which conditions the mechanism 
remains viable even when R parity is broken. 
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A deeper comprehension of the mechanism (origin, nature and scale) giving rise to the 
breaking of local supersymmetry (SUSY) still represents the major challenge we are facing 
on our way to construct realistic models of low-energy N=l supergravity models [|l|]. So 
far, the most conventional approach makes use of a so-called "hidden" sector |IJ which is 
made responsible for SUSY breaking at a large scale (10 10 — 10 11 GeV) with a gravitino 
mass of O(10 2 — 10 3 GeV) and the gravitational interactions representing the messenger 
of the SUSY breaking from the hidden to the observable sector. The alternative view that 
SUSY may be broken in a "secluded" sector at a much lower scale with gauge instead 
of gravitational forces responsible for conveying the breaking of SUSY to the observable 
sector had already been critically considered in the old days of the early constructions of 
SUSY models and has raised a renewed interest recently with the proposal by Refs. 0, |3|, [|] 
where some guidelines for the realization of low-energy SUSY breaking are provided. In 
these schemes, the gravitino mass (m 3 / 2 ) loses its role of fixing the typical size of soft 
breaking terms and we expect it to be much smaller than what we have in models with 
a hidden sector. Indeed, given the well-known relation |IJ between 777,3/2 and the scale 
of SUSY breaking yF, i.e. 7773/2 = 0(F/M), where M is the reduced Planck scale, we 
expect m 3 /2 in the keV range for a scale \/F~ of O(10 6 GeV) that has been proposed in 
models with low-energy SUSY breaking in a visible sector. 

In this letter we study some implications of SUSY models with a light gravitino (in the 
keV range) in relation with the dark matter (DM) problem. That such a light gravitino 
is very likely to be the lightest superparticle (LSP) and, hence, at least in models with R 
parity conservation, a good candidate for DM has already been known for quite some time 
now (since the early analysis of Pagels and Primack ||). The new point of our analysis is 
that there actually exist two different populations of gravitinos which are relic of the early 
Universe. First, we have gravitinos which through their 1/2 helicity component were in 
thermal equilibrium at some early epoch and, as we mentioned, are known || to be an 
interesting warm DM (WDM) candidate. Then, we have a kind of "secondary popula- 
tion" of gravitinos, which result from the decay of the next-to-the-lightest superparticle 
(NSP), presumably the lightest neutralino. They have a non-thermal distribution and 
exhibit features for the structure formation which are similar to those of a standard hot 
light neutrino in the tens of eV range. We will try to clarify in this letter under what 
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conditions these two populations of gravitinos can give rise to a tenable scheme of mixed 
DM (MDM). It will turn out that viable MDM realizations within the frame with light 
gravitinos that we envisage here lead to characteristic features both in the cosmological 
and particle physics contexts, making these models testable against astrophysical obser- 
vations and future accelerator experiments. In particular, on the astrophysical side, we 
find a relatively large He abundance (corresponding to slightly more than three neu- 
trino species), a suppression of high redshift galaxy formation with respect to the cold 
dark matter (CDM) scenario and a free-streaming scale of the non-thermal ("secondary") 
gravitinos independent of 7713/2, but sensitive to the NSP mass (with important conse- 
quences on the large scale structure formation). As for the particle physics implications, 
the implementation of a MDM scheme imposes severe constraints on the SUSY particle 
spectrum. For instance the lightest neutralino (NSP) should be an essentially pure gaug- 
ino and sfermions have to be rather heavy (in the TeV range). More specific features will 
be discussed below. 

As we previously mentioned, a stable particle with a mass in the keV range, like the 
light gravitino we are discussing here, was considered already long time ago as a possible 
warm DM candidate ||, i.e. a variant of hot DM but becoming non-relativistic at a 
much earlier epoch and, hence, having a much smaller free-streaming scale of 0(1 Mpc). 
Cosmological scenarios based on WDM were considered during the early '80s |J , soon after 
the shortcomings of HDM for large-scale structure formation were recognized. However, 
a purely WDM scenario seems to share the main drawbacks of standard CDM, once 
normalized on large scales to match the amplitude of the cosmic microwave background 
anistotropies J7J: (a) a spectrum of density fluctuations which is too steep on ~ 20 /i _1 Mpc 
f\ scales with respect to that measured for the distribution of galaxies and galaxy clusters 
J8|, and (b) too large fluctuations on ~ 10/i _1 Mpc scales, resulting in an overproduction 
of galaxy clusters 0. Colombi et al. have recently considered cosmological scenarios 
based on WDM. As a main conclusion, they found that a viable WDM candidate should 
have a mass-to-temperature ratio, m x /T x , twice that of the light neutrinos required by 



2 h is the Hubble parameter in units of 100kms _1 Mpc _1 ; 0.5^ h?z> 1 from observations; h = 0.5 is 
usually taken when considering critical density cosmological models with S7o — 1 , m order not to conflict 
with constraints on the age of the Universe. 
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the HDM model. For larger and smaller values of this ratio the WDM model rapidly falls 
into the CDM and HDM cases, respectively, thus requiring some degree of tuning for it 
to be a substantial improvement with respect of CDM. 

The difficulty for pure CDM frames has paved the way to a revival of the MDM 
|lT| where the cold and hot DM coexist, with a certain ratio of composites, in 



scenarios 



most cases of which the latter is supposed to be one (or more than one) light massive 
neutrino. In this scenario the presence of light neutrinos suppresses the growth of density 
fluctuations in the cold component on scales smaller than their free-streaming length. 
This goes in the right direction of generating a shallower spectrum, while keeping the 
fluctuation amplitude on the cluster mass scale to a more adequate level. 

However, it should be kept in mind that the cold+hot DM (CHDM) scenario is only 
one of the possibilities to implement the MDM idea. Another option that has been 



thoroughly investigated recently []T2j is that the hot thermal component is replaced by a 
volatile component made of particles with high rms velocity, which derive from the decay 
of a heavier particle. Such models based on cold+volatile DM (CVDM) differ from the 
more conventional CHDM schemes in that they involve a component which has a non- 
thermal phase space distribution function. An interesting implementation of the CVDM 
proposal is found in schemes where one and the same particle can play the twofold role of 
cold (or warm) and volatile component. An example was considered in Ref. ||13|| , where 
the fermionic partner of the axion, the axino, contributed the LSP. The decoupling of the 
axinos from the thermal bath yields the cold (or, rather, warm) component of the DM, 
whilst the axinos coming from the NSP decay constitute the volatile component which 
plays a role similar to that of light hot neutrinos as far as large scale structure formation 
is concerned. 

Here we point out that, thanks to a mechanism similar to that which was exploited in 
the abovementioned axino example, the light gravitinos can account for the warm as well 
as the volatile component, provided the NSP abundance before decay is large enough to 
yield a significant amount of non-thermal gravitinos. 

The helicity 3/2 component of gravitino has couplings of gravitational strength. If the 
Universe underwent the inflationary era, which we assume hereafter, its abundance was 
completely diluted during the inflation and is never produced later, so as to constitute a 
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significant portion of the mass density of the Universe. Thus the helicity 3/2 component 
plays no role in cosmology in this light gravitino case. On the other hand, the helicity 1/2 
component, or the longitudinal mode of the gravitino, has much stronger interaction [TJJ] 
when the gravitino is light having SUSY broken at a low energy scale. This is because 
it is essentially a Goldstino associated with the SUSY breaking, whose decay constant 
is given by the SUSY breaking scale \/~F ~ (r^^M) 1 / 2 . Indeed the explicit form of the 
interaction is [f0 



2 2 

£.// = a/T A tA ^ ^ xf ^ + ™if m tJ '**V + h - a > (1) 

where if) represents the helicity 1/2 component gravitino (the Goldstino) and m\, m x and 
are the masses of the gaugino, the chiral fermion and its superpartner, respectively. 
The helicity 1/2 component of the gravitino (hereafter we call it simply the gravitino) 
with such much stronger interaction than the gravitational one can be in thermal equilib- 
rium at an early epoch of the Universe. To see this, let us consider production/destruction 
of gravitinos a) by scattering (a + b «-> c + if)), and b) by decay and inverse-decay of a 
superparticle into a gravitino (a «-> b + if)). The total cross section for a + b — > c + if) was 



calculated in Ref. [0, being roughly 

htot -nml /2 M^ {Z) 
where m- g is the gluino mass. The resulting interaction rate is 

r 'scatt — ^tot n rad (3) 

with n ra d — C(3)/7r 2 7" 3 being the number density for one massless degree of freedom. 
Comparing this with the expansion rate of the Universe H, one finds that the interaction 
is effective to keep the equilibrium (T scatt ^ H) until the temperature becomes^ 



(4) 



3 For T?S nig the above estimate for S tot will not be accurate. However, it will not change our argument 
drastically. 
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Next consider the decay and inverse-decay process. The decay width of a superparticle 
(R-odd particle) a into its superpartner b and a gravitino is given by |L7|, [L6| 

where m a is the mass of a. Suppose that, for simplicity, all superparticles in the MSSM 
have the same mass, iris- Given a temperature T > ms, the interaction rate for this 
process is estimated as 

r decay = (T) x 

~ g s —Y[a -> &V) 
487rm| /2 M 2 T' 



(6) 



Here represents the effective degrees of freedom of the superparticles in the MSSM, 
~ 100, and ns is the number density of the superparticles. From this it follows that 

^decay _ ffg m S f»\ 

H ~ 50^(T)V2 m 2 /2MT 3' ^ > 

where g*{T) represents the effective massless degrees of freedom. Thus one may conclude 
that for O(10 2 )GeV(m3/2/lkeV) 3//2 , the decay and inverse-decay process is effective 
to maintain the equilibrium of the gravitinos as long as the superparticles are relativistic. 
As they become non-relativistic, the interaction rate drops exponentially, ~ e _ms//T , due 
to the Boltzmann suppression of the number density of the superparticles. In order to 
determine the freeze-out temperature Tf precisely, one must explicitly take into account 
the superparticle mass spectrum when integrating numerically the Boltzmann equation. 
Here we will not attempt to do so, rather we just estimate that the effective massless 
degrees of freedom at Tf is g*(Tf) ~ 100 — 200, admitting a factor 2 ambiguity. This is 
enough for the purpose of this paper. 

Suppose that the gravitinos were once in thermal equilibrium and were frozen out 
with g*{Tf) given above. We assume that the gravitino is the LSP and stable, which is 
guaranteed by R parity conservation. We will discuss later the case where R parity is not 
conserved. Following the standard procedure, the density parameter Q t h contributed by 
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relic thermal gravitinos is 



n t ,h 2 = 1.17 



( m 3/2 \ ( g*(T f Y 
UkeVy I 100 



-i 



or 



m 3/2 = 0.85keV (n th h 2 ) {^^j ■ (9) 

Therefore, a gravitino in the abovementioned keV range provides a significant portion 
of the mass density of the present Universe and behaves as a WDM candidate, as was 
mentioned previously. 

We now turn to another contribution of the gravitino abundances, namely gravitinos 
produced from massive particle decays after the freeze-out. To be specific, we consider the 
gravitinos produced by the non-thermal decays of binos, the superpartner of the U(1)y 
gauge boson, assuming that it is the NSP, i.e. the lightest among the superparticles in 
the MSSM sector. 

^From Eq. (f|), one finds that the decay width of a bino into a photon and a gravitino 

is 



T(6->7^) = cos 2 9 w x 



1 m l 
48tt m\ /2 M 2 



2.15 x lO^GeV f-^V f^Y 



(10) 



V30GeV/ \m z/2 ) ' 

where m~ b is the bino mass and 9w stands for the weak mixing angle. 

Equating the above decay rate with the expansion rate of the Universe, one determines 
the bino decay temperature T D : 



(11) 



Hence the bino decay takes place well before nucleosynthesis starts. 

Each bino produces one gravitino at decay. From this, it follows that the abundance 
of the non-thermal gravitinos coming from the bino decays is 



O _ m 3/2 , , 

* L non-th — iL b K 1 ^) 
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where represents the would-be bino mass density at present, if it were stable. For 
example, for 777.3/2 = 0.5 keV, m~ b = 30 GeV, one needs Q~ b = 1.2 x 10 7 to achieve fl non - th = 
0.2. 

This large abundance required for the next-to-the-lightest superparticle (NSP) imposes 
severe constraints on mass spectrum of superparticles. First of all, the NSP must be 
almost pure gaugino: a non-negligible contamination of higgsino component would allow 
it to annihilate through, for example, s-channel Z boson exchange so that one would not 
expect such a large f^. Second, it requires that sfermions (i.e. squarks and sleptons) must 
be very heavy. The relic abundance of the bino is estimated as 

(m| + m~) 2 

ttzh 2 ~ 1 x 1(T 6 f 9 b GeV~ 2 , (13) 

771- 
b 

where we have assumed for simplicity that all sfermions have the same mass mj. For 
m~ b = 30 GeV and h = 0.5, the sfermion mass mj must lie around 7 TeV to obtain the 
abovementioned value Q b = 1.2 x 10 . 

This mass pattern provides a non trivial constraint for a model of low-energy SUSY 
breaking. In the model of Ref . , the messenger sector connects the visible sector with 
the supercolor sector in such a way that gauginos can acquire their masses at one-loop 
level whereas squarks and sleptons get mass squared at two-loop level. The ratio between 
the gaugino and sfermion masses is model-dependent. Indeed, it is strictly related to 
the specific pattern of the mass matrix of the messenger fields. One can envisage so- 
lutions where a mass pattern suitable for the scheme that we propose here is actually 
implemented.^ 

Thus far, we have obtained two different populations of the gravitino DM. One is the 
thermal gravitino which gives the WDM, and the other is the non-thermal gravitino whose 
properties are similar to the HDM. To clarify the latter point, we calculate the redshift 
Z nr at which the non-thermal gravitino becomes non-relativistic. The momentum of the 
non-thermal gravitino gets red-shifted as 

WW = ^ (w) 



J We thank G. Giudice for discussions on this point. 
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where a(t) is the expansion factor at a given time t and is its value at the bino decay. 
Let a nr be the expansion factor when the gravitino becomes non-relativistic, i.e. when 
kphys = 772.3/2- From Eq. (|14D, it follows that 

m- b a D 

m 3/2 = — • (15) 

_ (X 



nr 



Therefore, it turns out that 

Znr + 1 



CIq CIq Q,£) 

Uj nr 

1/3 



9*s(Tp) \ T D 2m 3/2 
a*s(To) I T m ; 



b 



' g.(TD) Y' 12 ( m- b Y' 2 



= 6-42 x 10^^ j ^ . (16) 

For a heavier m^, the temperature at its decay Tjj will be larger than the freeze-out tem- 
perature of the bino, ~ m~ b /20 and then the above argument for Z nr should be modified. 
We will not detail this here, since a heavier yields a larger Z nr , which is presumably 
less attractive from the viewpoint of large scale structure formation. 

It is interesting to compare this with the redshift at which the thermal gravitinos get 
non relativistic. Since Z is related to the gravitino temperature T 3 / 2 according to 

z | -i _ fl o _ ( 9*{T f ) \ 1/3 T 3/2 
Z + l -^-\^m) T ' 

and the thermal gravitinos become non-relativistic around when their temperature be- 
comes 7713/2/3, it turns out that 

Znr . „ I 9*(T f ) \ 1/3 m 3/2 /3 
\9*s(T ) J T 

= 4.14 x 10 6 x ( 9 -^A) 1/3 (18) 
V 100 J VlkeV/ V ; 

Hence thermal gravitinos become non-relativistic much earlier than the non-thermal 
ones. 

Once Z nr is known, one can estimate the free streaming length until the epoch of the 
matter-radiation equality, Xfs, which represents a quantity of crucial relevance for the 
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formation of large-scale cosmic structures. If v(t) is the typical velocity of a DM particle 
at the time t, then 



a(t) 

= 2t x -^[l + \ n {Z nr /Z eq )} 

= 6.08 x 10 5 x Z^Mpcfl + ln(Z nr /2.32/i 2 x 10 4 )] (19) 

According to this estimate, for the non-thermal gravitino with Z nr = 6.42 x 10 4 , the 
free-streaming length is A ~ 30 Mpc, which corresponds to a supercluster size. Thus the 
non-thermal gravitino in our scenario will exhibit properties similar to those of a regular 
HDM candidate, like a light thermal neutrino, as far as large scale structure formation 
is concerned. On the other hand, the free-streaming length for the thermal gravitinos is 
about IMpc (for Z nr ~ 4 x 10 6 ), which in turn corresponds to ~ 1O 12 M . 

Therefore, this scenario corresponds to a MDM model, with warm and volatile DM 
components provided by thermal and non-thermal gravitinos, respectively. A similar 



scheme has been recently considered by Malaney et al. [H5|, where the hot component 
is represented by ordinary light thermal neutrinos, while the warm part is provided by 
sterile neutrinos. These authors conclude from their analysis that this model is virtu- 
ally indistinguishable from the cold+hot DM scenario as far as the large-scale structure 
formation is concerned. A crucial difference between this scheme and the one we are 
proposing here lies in the fact that for the thermal neutrinos the value of Z nr is fixed by 
their mass and, therefore, by the contributed density parameter Q u . On the contrary, Z nr 
for the non-thermal gravitinos does not have a one to one correspondence with fl n on-th 
[cf. Eq. QT§)]. 

As for the warm component, taking 777.3/2 ~ 200 eV and being 

T:i/2 K (i?^y\ (20) 

\9*{ T d)J 

the relation between the neutrino and the gravitino temperature, it turns out that the ratio 
777.3/2/T3/2 is at least 10 times larger than m v /T V) thus showing that the warm component 
behaves like CDM, at least on scales ^ 1 /i _1 Mpc. 
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Pierpaoli et al. |12] have recently considered a MDM scheme, in which thermal axinos 
play the role of CDM, while non-thermal axinos provide the volatile part. Therefore, as 
far as the cosmological implications are concerned, this model behaves like the one based 
on gravitinos, at least on mass scales larger than that of a galaxy. As a main result of their 
analysis, Pierpaoli et al. [|12| pointed out that the simultaneous request of reproducing 
the observed abundance of galaxy clusters |§ and of high-redshift (Z ~ 4) damped Ly-ct 
systems (DLAS) |20| requires Q non - th ~ 0.2 and Z nr ^, 10 4 . As for the cluster abundance, 
since it is determined by the fluctuation amplitude on scales ~ 10 /i _1 Mpc, we expect 
similar predictions when replacing the cold component with the warm one. However, since 
DLAS are ought to be associated with protogalaxies, the constraint they provide refers to 
scales ^ 1 /i _1 Mpc. In this regime, the effect of free-streaming in the warm component 
suppresses the fluctuation growth, therefore decreasing the predicted abundance of high- 
redshift DLAS. This is potentially a critical test for our DM model, the DLAS abundance 
being already recognized as a non-trivial constraint for usual cold+hot DM models f2~lf . 
It is however clear that more quantitative conclusions would at least require the explicit 
computation of the fluctuation power spectrum for our warm+volatile DM model, which 
is beyond the scope of this letter. Furthermore, before definitely assessing the confidence 
level at which such observational constraints rule out a model, a further clarification of 
both the reliability of available data and of the corresponding interpretative framework 
are required. 

As already pointed out, the value of Z nr for non-thermal gravitinos in principle does 
not depend on f2 non _^. On the other hand, a large amount of exotic relativistic particles 
around the era of the big-bang nucleosynthesis would contribute to the energy density of 
the Universe at that epoch, so as to accelerate the expansion of the Universe and result in 
a significant increase of the 4 He abundance. In our case, both thermal and non-thermal 
gravitinos will contribute. It is convenient to express such contributions in terms of the 
effective number of extra generations of neutrinos defined by AN U = Ap/p u , where Ap is 
a contribution to the energy density by an exotic particle and p v is the energy density of 
a neutrino (one species). The contribution from the thermal gravitino is easily obtained 
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as 



On the other hand, to evaluate the contribution from the non-thermal gravitino, one 
should note that the energy density of the gravitino p non -th evolves differently from that 
of the neutrino p v only after the gravitino becomes non-relativistic: 

Pnon-thi^o) Q"nr 



^pjnon-th _ Pnon-th 
Pu 

^From this it follows that 



(22) 

lMeV Pv{to) a 



AJC~* = 0.134 x (°~* > h2 ) ( 6 42 X W ') . (23) 
\ 0.05 J \ Z nr J 

The total contribution of the gravitinos to the energy density during the nucleosynthesis 
is the sum of Eqs. ( pT|) and fl23|). Note that it is comparable to, or even larger than the 
2a upper bound for AN U = 0.16 coming from the observations of the 4 He abundances 
obtained recently by Olive and Scully p2[ . As was cautioned by these authors, the sys- 
tematic error they used to obtain this number might be somewhat underestimated. On 
the other hand, Kernan and Sarkar claimed that up to AN V ~ 1.5 is allowed by 
observations of high D abundance in high-redshift Lya clouds, while no general agree- 
ment exists about deuterium observations in high-redshift QSO's (see, e.g., ref . p4 and 
references therein). Keeping this in mind, we believe that it is premature to exclude our 
scenario by this constraint of the 4 He abundance. Rather we should emphasize that our 
scenario requires a significant excess of the primordial 4 He and D abundance which might 
be excluded or confirmed in future. 

So far, we have assumed that R parity is strictly conserved. As we will see soon, the 
light gravitino is long lived enough to account for DM even when R parity is broken. 
This contrasts with the case where the LSP is a neutralino: unless R is violated by an 
extremely tiny amount, one has to demand R conservation to ensure that the LSP is a 
viable DM candidate. 

R parity breaking inevitably leads to baryon or lepton number non-conservation. 
Given the lightness of the gravitino that we consider here, it cannot decay into baryons. 
Hence it is stable if R is violated in the baryonic sector. Next we will estimate the lifetime 
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of the gravitino when the lepton number is not conserved. To be specific, we concentrate 
on the case that a term Xij^LiLjE^. with a Yukawa coupling appears in the super- 
potential as an R violating interaction, where L and E c are chiral multiplets of SU(2)l 
doublet and singlet leptons, respectively. The Latin subscripts represent generations and 
no summation over them is taken. At the one-loop, the gravitino decays to a photon and 
a neutrino (or an anti-neutrino) through lepton-slepton loops, with the decay rate 

r(V> -> u(i?)j) = ^a^ 2 [K (ln(m fi ./m/j) 2 + m\ (\n{m t Jmi k )f 

+ mf. (\n{m~ lR Jm h )f +mj (\Yi{m~ lR Jm h )f , (24) 

where a is the fine-structure constant, a>\ = Xfj k /4:Tr, mi the mass of the charged lepton, 
and mj L (mj- ) is the mass of the left-handed (right-handed) slepton. The decay rate will 
be maximized when Aj33 7^ (i = 1, or 2). If we assume for simplicity that all charged 
sleptons have the same mass mr = mr = mj, then the decay rate reads 

r(^ -> z/(z/)7) = 32?r 3 (ln(m r /m r )) (25) 

For mj = 7 TeV, one finds the lifetime 

T y, _> „(p) 7 ) = 2 x lO 21 ^ 1 (^) _1 sec, (26) 

which is larger than the present age of the Universe. 

The radiative decay of such a long-lived DM would contribute the diffuse extragalactic 
photon background ||26|| . The photon number flux emitted by the gravitino is estimated 



as 

for ft, = 0.5 and we have assumed that the total gravitino mass density saturates the 
critical density of the Universe. Observations of the diffuse photon background put a 
constraint on the lifetime 

r> 10 25 - 10 26 sec (28) 

for m 3 / 2 — 0.3 — IkeV, implying «a~ 10~ 5 . Therefore the gravitino with lifetime longer 
than Eq. (ESf) will be accounted for a viable DM. 
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The abundance of the non-thermal population of the gravitino may be affected by the 
R parity violation. A necessary condition to avoid this failure is that the NSP (the bino 
in our case) dominantly decays to the gravitino, not to R even particles. To illustrate how 
this condition puts constraints on R parity breaking interactions, consider again the term 
XijkLiLjEl in the superpotential. Through this interaction, the bino can decay to three 
leptons, one (anti-)neutrino and two charged leptons. We can estimate this R violating 
decay rate as 

x 4- (29) 



aa x m~ 

V — 

1927T cos 2 9w mf 



i 

We then require that Eq. (p9| ) should be smaller than the partial decay width to a gravitino 
and a photon given by Eq. ([TUP, yielding a constraint 

4 

Tfl~ 

«a<4cos 4 ^-^— , (30) 

which is 0(1CT 7 ) for the parameter range we are considering. 

To conclude, we have considered the case where the light gravitinos frozen out from 
the thermal bath constitute a warm DM, whereas the non-thermal gravitinos produced 
by the NSP decay contribute as a volatile component of the DM. This mixed DM sce- 
nario requires a certain pattern of superparticle mass spectrum, which may be tested in 
future collider experiments. From the astrophysical side, the main implications of our 
DM scenario can be summarized in the two following points: (a) as for the primordial 
nucleosynthesis, the effective extra neutrino generations associated to gravitinos imply a 
rather large abundance of 4 He and D; (b) as for large-scale structure formation, the resid- 
ual free-streaming of the "warm" thermal gravitinos delay the galaxy formation epoch 
with respect to the "cold" case, so that the observed high-redshift structures are required 
to be associated to dwarf protogalaxies. 
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